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Scanning tunneling microscopy (STM) measurements have observed modulations in the density 
of states (DOS) of a number of high-T c cuprates. These modulations have been interpreted in 
terms of electron-boson coupling analogous to the dispersion "kinks" observed by angle-resolved 
photoemission spectroscopy (ARPES). However, a direct a reconciliation of the energy scales and 
features observed by the two probes is presently lacking. In this paper we examine the general 
features of el-boson coupling in a d— wave superconductor using Eliashberg theory, focusing on the 
structure of the modulations and the role of self energy contributions X z and A^. We identify the 
features in the DOS that correspond to the gap-shifted bosonic mode energies and discuss how the 
structure of the modulations provides information about an underlying pairing mechanism and the 
pairing nature of the boson. We argue that the scenario most consistent with the STM data is that of 
a low-energy boson mode renormalizing over a second dominant pairing interaction and we identify 
this low-energy mode as the out-of-phase bond buckling oxygen phonon. The influence of inelastic 
damping on the phonon-modulated DOS is also examined for the case of Bi2Si2CaCu208+i . Using 
this simplified framework we are able to account for the observed isotope shift and anti-correlation 
between the local gap and mode energies. Combined, this work provides a direct reconciliation of 
the bandstructure renormalizations observed by both ARPES and STM in terms of coupling to 
optical oxygen phonons. 

PACS numbers: 74.25.Kc, 73.40.Gk, 74.72.-h 



I. INTRODUCTION 

In conventional superconductors, the role of electron- 
phonon (el-ph) interactions as the pairing "glue" was 
confirmed by the observation of fine-structure in the 
tunneling-derived density of states (DOS). 1,2 In the 
cuprates, a great deal of effort has been expended in 
search of similar signatures using different probes. 3-35 
The observation of a particular mode, be it spin or lattice, 
might be used to infer the nature of the principal pairing 
mediator. However, to date, there is no consensus on the 
existence of such a mode 36 and if it does exists, whether 
it is tied to the lattice or spin degrees of freedom. 

Electronic renormalizations in the form of dispersion 
"kinks" have also been observed in the electronic dis- 
persion of the cuprates by angle-resolved photoemission 
spectroscopy (ARPES). 6 ~ 8 ' n < 13 - 18 Occurring at an en- 
ergy scale ~ 65 — 70 meV in the nodal direction (0,0) - 
(7r,7r) of the Brillouin zone, these kinks have been in- 
terpreted as coupling to a bosonic mode of either an 
electronic 6-12 or lattice origin. 13-21 According to Eliash- 
berg theory, the position of a kink due to coupling to a 
sharp bosonic mode in a c£-wave superconductor is ex- 
pected to occur at an energy of Ao + fioj where Q,q is 
the energy of the mode and Ao is the maximum value of 
the d-wave gap 2,16. With this observation, the energy of 
the mode responsible for the kink has been estimated at 
^35 — 40 meV. This energy coincides with the energies 



of the spin-resonance mode centered at Q = (n/a,iT/a) 
observed by neutron scattering 37 and the out-of-plane 
"Big" bond buckling oxygen modes. 15 As a result, there 
is considerable debate as to the identity of the underlying 
mode, however the observation of multiple mode coupling 
in some cuprates, 20,21,28 as well as recent measurements 
of an isotope effect for the kink energy, 14 provide strong 
evidence for the phonon interpretation of the ARPES 
renormalizations. However, regardless of the identity of 
the boson, the ~ 35 — 40 meV energy scale extracted from 
ARPES data appears to be in conflict with the ~ 52 meV 
energy scale extracted from STM by Lee et al.. 3 There- 
fore it is an open question as to whether the two experi- 
mental probes are observing different mode couplings or 
if they are reflecting different manifestations of the same 
bosonic mode. 

Microscopic inhomogeneities in the local density of 
states (LDOS) which is proportional to the derivative 
of the tunneling current N(u) oc dl/dV, as well as signa- 
tures of coupling to a bosonic mode, have been observed 
in scanning tunneling microscopy (STM) experiments on 
Bi 2 Sr2CaCu20 8+(5 (Bi-2212). 3 - 5 In Ref. 3, estimates for 
the local gap size are determined from the peak-to-peak 
distance of the coherence peaks while the energy of the 
bosonic mode is identified as the energy position of a peak 
in d 2 I/dV 2 measured relative to the energy of the super- 
conducting gap. While the positions of the supercon- 
ducting coherence peaks vary at different tip locations, 
estimates for the mode energy are inversely correlated 
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with the local gap size and the distribution of mode esti- 
mates, centered at ~ 52 meV, shows a clear isotope shift 
upon 18 O substitution. The energy of the bosonic mode 
also appears to be immune to doping while the spectra 
changes qualitatively. These observations are inconsis- 
tent with a coupling to the spin resonance mode, 35 and 
points to a lattice origin for the mode involving oxygen 
vibrations. 

Fine structures in the DOS of a number of additional 
cuprates have been reported by other STM and SIS 
junction tunneling measurements, each producing dif- 
ferent estimates for mode energies. A structure simi- 
lar to that reported for Bi-2212 has been reported in 
Bi 2 Sr2Ca2Cu30io +(5 (Bi-22 23), 22 with the energy scale 
of the mode estimated at ~ 35 meV. It is important to 
note that in Ref. 22 the mode energy was associated 
with the position of the minima in the LDOS (root in 
d 2 I/dV 2 ) relative to the local gap size as opposed to 
Ref. 3, which extracted the mode estimate from the 
peak in dI 2 /dV 2 . Bicrystal grain boundary SIS junction 
measurements on optimal doped Lai.84Sro.i6Cu04 thin 
films observe modulations in the DOS which correspond 
well with peaks in the neutron derived phonon spectra, 23 
and in agreement with the multiple features present 
in ARPES data. 28 Examinations of tunneling data on 
YBa2Cu307_(5 (YBCO) and the electron doped system 
Pro.88LaCeo.i2Cu04 have produced a similar correspon- 
dence between modulations in the DOS and phonon den- 
sity of states in these materials. 24 ' 25 The observation of 
multiple mode coupling, as well as the fact that the spin 
resonance is well separated in energy from the phonons 
in the electron doped systems provide further evidence 
that the spin resonance mode is unlikely to be the source 
of these features as proposed by Ref. 35. However, the 
question remains whether coupling to the mode can be 
cast in the usual form for el-ph coupling to oxygen modes, 
as inferred from ARPES 15 ' 19 , or whether the structure in 
the LDOS could be due to phonon-assisted co-tunneling 
from the tip via the apical atom. 27 

Co-tunneling via a strong local coupling of electrons in 
the STM tip to the apical atom imparts structure in the 
form of peaks in the dl/dV at an energy scale of A + fi 
and multiples of the phonon frequency Ao + 20, Ao + 30, 
... even though the coupling of the planar superconduct- 
ing electrons to the apical atom may be weak. 27,38 How- 
ever, recent observations of charge-ordering structures 
via Fourier transform spectroscopy in Bi-2212, 39 which 
has apical oxygen atoms, and Ca2_ a; Na 2: Ca02Cl2, 40 
which doesn't, empirically indicate that no one particu- 
lar matrix element controls the tunneling pathway along 
the c-axis. If a single pathway were to dominate the 
tunneling process one would expect the charge ordering 
features in these two materials to have qualitatively dif- 
ferent structures. 

In light of these observations we investigate the role 
of el-ph coupling between planar oxygen vibrations and 
the electrons of the Cu02 plane. The aim of this paper 
is twofold. First, due to the discrepancies in the struc- 



tures in dl/dV used to extract boson mode energies in 
previous works, we examine the qualitative signatures of 
el-boson coupling in a <i-wave superconductor. In doing 
so, we explicitly determining which features in dl/dV (or 
d 2 I/dV 2 ) are best identified with the energy of the boson 
mode. Furthermore, we demonstrate how the qualitative 
structure in dl/dV can be used to determine if the mode 
coupling is the predominant pairing mechanism. We ar- 
gue that most consistent interpretation of the STM data 
is that of a low-energy mode renormalizing over a domi- 
nant high-energy (or instantaneous) pairing interaction. 
In this case, the mode energy should be extracted from 
a minimum in N(u>) cx dl/dV. In light of this find- 
ing, the mode energy estimate from Ref. 3 is revised 
to ~ 35 — 45 meV, in agreement with Ref. 22 and the 
energy scale extracted from ARPES measurements. 16 ' 21 
The second goal of this paper is to demonstrate that the 
modulations in N(w) can be understood in terms of the 
same el-ph coupling models that have been successful in 
accounting for the dispersion "kinks" observed ubiqui- 
tously in the cuprates by ARPES. Combined, these cal- 
culations demonstrate a unified picture of el-ph coupling 
in the cuprates probed by both STM and ARPES. 

The organization of this paper is as follows. In sec- 
tion II the Eliashberg framework for calculating the self- 
energy due to el-boson coupling in a cZ-wave superconduc- 
tor is reviewed. In section III focus is placed on el-boson 
contributions in multiple momentum channels in order 
to demonstrate the qualitative changes in the self-energy 
that are expected for a mode which couples differently in 
each of the momentum channels. We find that these con- 
siderations are critical and can qualitatively change the 
structure of the boson modulations in the DOS, alter- 
ing the feature in dl/dV corresponding to the mode en- 
ergy. This aspect of the calculation has generally been ne- 
glected in previous Eliashberg treatments, 12,29-31 where 
the boson contribution to the single-particle (A z ) and 
anomalous (A^) self-energies have been taken to be equal 
to (or proportional to) one another over the entire en- 
ergy range of the boson spectrum. For an s-wave su- 
perconductor X z = A<£ but for a d-wave superconductor 
\ z ^ \<j, and usually X z 3> A^ in modes derived from 
the strong Coulomb repulsion having strong onsite (X z ) 
and near-neighbor (A^) interactions. Next, having estab- 
lished the role of the symmetry channels, we then turn 
to the qualitative differences in the el-boson structures 
in the infinite and finite band formalisms of Eliashberg 
theory. Since the cuprates have a relatively narrow band- 
width due to the Cu 3d character of the pd-a* band, 
one expects that the latter case is the more appropri- 
ate formalism for these materials. Finally, in section IV 
we present a model calculation for the el-ph modulated 
LDOS in Bi-2212. Here it is demonstrated that the DOS 
modulations can be reproduced using coupling strengths 
similar to those used in previous ARPES treatments. 15 
We then introduce a phenomenological treatment of local 
damping effects reflecting the inhomogeneity observed in 
Bi-2212 and the associated broadened spectral features. 4 
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This model is able to reproduce the anticorrelation be- 
tween the local superconducting gap estimates and local 
mode estimates. Finally, the reported isotope shift of the 
distribution mode estimates is naturally captured by the 
el-ph coupling model, as our calculation explicitly show. 
We then end with a brief summary and discussion about 
the implications of this work. 



II. ELECTRON-BOSON COUPLING 

In Migdal-Eliashberg theory, which neglects crossing 
diagrams, the self-energy due to el-boson coupling is 
given by 



1 f ' 

Z(k,iu n ) = —J2Jdis- 



2va 2 F(k,( l , v) 



q,m 



/ 2 + {uj m - 0J n ) 2 

f 3 G(p, iu! m )T 3 (1) 

where w m (w„) is a Boson (Fermion) Matsubara fre- 
quency, and a 2 F(k, q, v) = — \g(k, q)| 2 ImZ?(q, v) is the 
effective electron-boson spectral function. Here, the ver- 
tex g(k, q) gives the strength of scattering of the electron 
from state k to p = k q and D(q, icu m ) is the boson 
propagator. 

In the superconducting state, the electron propagator 
G(k, iuj n ) is given by 

G -1 (k, iuj n ) = iuj n Z(k, iuj n )r + [e k + x(k, iw„)]f 3 

+ 0(k,iw„)fi. (2) 

Here, the self-energy has been divided into the canonical 
form 2 

£(k,iw„) = iuj n [l - Z(k, ico n )]to + x(k,iw n )^3 

+ 4>(k, iw n )Ti (3) 

where iuo n [l — Z(k,iuj n )\ and x(k, ioj n ) are the odd and 
even components of the electron self-energy, respectively, 
and 4>(k, iuj n ) is the anomalous self-energy. Eqs. (1) - 
(3) can now be self-consistcntly solved on the Matsubara 
frequency axis. An analytic expression for the self-energy 
on the real axis can be derived from these equations by 
explicitly performing the analytic continuation iuj n — > 
lu + i8, where 5 > is an infinitesimal real number. The 
single-particle DOS is N(co) is obtained from the electron 
Green's function N(w) = — 2^ k ImGn(k,a;)/7r. 

Marsiglio et al. 41 have developed an efficient iterative 
procedure for performing the analytic continuation of the 
el-boson self-energy. In this formalism S(k, uj) is ob- 
tained by iteratively solving 41,42 (p = k q) 

S(k,w) = — A o(k,q,w - icj m )f 3 G(p,icj m )f 3 



q m— 



[n b {v) + n f (v - w)]f 3 G(p, lo - v)t 3 



(4) 



where rib and rif are the Bose and Fermi factors, respec- 
tively, and 



A (k, q,w) = / dva 2 F(k, q, v)- 



2v 



In evaluating Eq. (4), the solution to the imaginary 
axis equations S(k, iuj n ) are used as input and S(k,cj) 
is solved for iteratively. 41 

In large bandwidth systems the momentum sum is typ- 
ically evaluated by assuming a linear band near the Fermi 
level £k = (k — hp) ■ vp (with vf the Fermi velocity), re- 
placing the normal state density of states N(uj) with its 
value at the Fermi level ./V(O) and extending the limits 
of the energy integral to infinity. 41 (Note that in this ap- 
proximation x(k, uj) is featureless and only contributes to 
an overall shift of the chemical potential. 2 ) For brevity 
we refer to this approximation as the "infinite band" for- 
malism of Eliashberg theory. 



III. QUALITATIVE SIGNATURES OF 
EL-BOSON COUPLING 

Before turning to the model calculations, let us dis- 
cuss the qualitative signatures of el-boson coupling in a 
(i-wave superconductor. We would like to examine cou- 
pling in various momentum channels and, to this end, the 
coupling constant can be expanded in terms of Brillouin 
zone harmonics 



Igik,^ 2 =J2Yj(k)\g J ,j,\ 2 Y J ,(k-< l ) 



(5) 



J, J' 



where the sum J runs over the irreducible representa- 
tions of the point group of the crystal. We assume gj^ji 
is diagonal in this basis and identify J — and J = 2 



with the s and cL 



symmetries. If one then only ad- 



mits the gap solution in the J = 2 channel, the energy 
and momentum dependence of the self-energies are fac- 
torable with Z{k,u) = Z(uj)Y {k), x(k,w) = X {u)Y (k) 
and (f>(k, uj) — (f>(uj)Yj (k) , where J = 0,2 for each sym- 
metry. 

The overall strength of the el-boson coupling in mo- 
mentum channel J can be parameterized by the dimen- 
sionless constant Aj 



A, = 



2du E k , q a 2 F (k, q, iy)Yj(k)Yj(p)S(e ii )S(e p ) 
E k ^(k) 2 <5(e k ) 



(6) 

It is important to note that A,/=o = A z {Yo(k) = 1) char- 
acterizes the contribution to the el-boson self-energies 
Z(k,co) and x(k, ui) while \j = 2 = ^<p (Y 2 = [cos(k x a) — 
cos(k y a)]/2) characterizes the contribution to <j)(k,u) for 
a d x 2_ y 2 superconductor. The relative values of \ z ,<p a l so 
determine the transition temperature T c , which, in the 
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FIG. 1: (Color online) (al), (a2) The DOS N(w) and corre- 
sponding real and imaginary parts of the gap function A(w), 
respectively, for an s-wave superconductor coupled to a sin- 
gle bosonic mode, (bl), (b2) N(ui) and A(uj) for a d-wave 
superconductor with A z = 1.6 and A^, = 0.8A Z . The insets 
of panels (al),(bl), show dN/duJ and a 2 F(y + Ao) in order 
to highlight the correspondence between features in the DOS 
and peaks in the boson spectrum. The secondary structure 
in A(lu) at uj = Ao + 2Q.Q are due to self-energy effects due to 
two-phonon processes. 



limit of weak coupling is given by 



kbT c = 1.13M7o ex P 



A, 



(7) 



The relative magnitudes of \ z ,4> can also affect the quali- 
tative signatures of the boson modulations in the density 
of states. 

In the superconducting state, neglecting the role band- 
structure, the DOS N s (uj) can be written as 



N n (0) 



Re< 



\y/0J 2 - A 2 (k,w) 



(8) 



where N n (0) is the density of states at the Fermi level 
and A(k,w) = 0(k, u)/Z(k, uj) = Ai(k,w) + iA 2 (k,u;) 
is the complex momentum-dependent gap function and 
(...) denotes an average over the Fermi surface. 2 In what 
follows, the notation Ao is introduced for the maximum 
value of the superconducting gap on the Fermi surface, 
which is defined by the value of the gap function at the 
gap edge Ao = A(k^ N ,w = Ao), where k^ N is the Fermi 
momentum along the zone face (0,0)-(tt,0). 

For uj ^> Aq, Eq. (8) can be expanded yielding 2 



N s (uj) 
N n (0) 



= l + ^(A 2 (k, w )-A 2 (k, w )). (9) 



From Eq. (9) it is clear that the phononic substructure is 
given by the frequency dependence of the A(k, uj), which 
is obtained by evaluating Eq. 4. In Figs. Ial-a2 we plot 



N s (ui) and the corresponding A(ui), respectively, for an 
isotropic s-wave superconductor. Here, we have assumed 
a single bosonic mode with a Lorentzian spectral density, 
characterized by a half- width at half-maximum of Ti, = 1 
meV and centered at ilo = 52 meV (inset of panel Fig. 
lal). The overall coupling strength has been set to A z = 
0.8. 

As uj — > Ao + f^o, the real part of A(w) begins to rise 
producing an enhancement in the DOS for uj < A + 
H,q. At A + r^o the real part begins to drop while the 
imaginary part experiences a sudden rise due to Kramcrs- 
Kronig consistency. This results in a rapid suppression 
in the DOS at this energy scale, which drives the DOS 
below its bare value. As a result, the energy scale Ao + £!o 
manifests as a shoulder in N(uj) or a minimum in dN/duj, 
as shown in the inset of Fig. lal. This is the classic 
McMillan-Rowell signature, similar to that observed in 
Pb, where phonons are solely responsible for pairing. 1,2 
In the case of a d-wave superconductor the situation is 
nearly identical when A z = \<f, = 0.8, and the magnitude 
of the gap function is comparable to that obtained for 
the s-wave case with a similar value of A z (not shown). 
Again, since the contribution from the boson is equal 
in the two channels, the boson energy scale manifests 
as a shoulder on the low-energy side of the modulation. 
This signature qualitatively unchanged when A z > A^, 
as shown in Figs. Ibl,b2 for A z = 1.6, A^ = 0.8A Z , with 
the overall magnitude of A(ui) reduced by a factor of 2. 
The reduction in A(u>) illustrates the importance of the 
relative magnitudes of A z and A^. 

We now consider the case of a low-energy boson which 
renormalizes over a second mode, higher in energy and 
dominant in its contribution to pairing. We associate this 
high energy mode with spin fluctuations with a broad 
Lorentzian spectral density centered at Q s f ~ 2 J = 260 
meV and with T s f — 30 meV. The total coupling to 
this mode is taken to be \ s f,z = ^sf,<t> — 1-6. For the 
low-energy boson we take Ft = 0.5 meV, A z = 1.6 and 
A = O. 

The results for the two-mode calculation are shown 
in Fig. 2. The structure of the renormalizations in 
this case differs considerably and the shoulder on the 
low energy side of the renormalization is significantly 
less pronounced. Here, the sharper spectral density and 
increased value of A z for the low-energy mode are re- 
quired in order to accentuate the weak low-energy fea- 
ture. Without this increased coupling the low-energy 
renormalizations are difficult to resolve. In the Eliash- 
berg formalism, a high energy boson produces a gap 
function whose real part is relatively frequency indepen- 
dent for energies on the order of lower energy mode f2{,, 
while an instantaneous pairing interaction produces a fre- 
quency independent gap function up to an energy scale 
set by the Coulomb interaction. 43 Therefore, assuming 
a frequency independent pair field <po, modulated by el- 
boson coupling, we write 

(j)(u>) + 6<f>(ui) (j> 1 + 



Z(ui) Z + 5Z{ui) Z a 1 + 5Z/Z 
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FIG. 2: (Color online) (a) N(lj) for a d-wave superconductor 
where the £7 ~ 52 meV mode {Xt,z = 1-6, Xb,4, = 0) renormal- 
izes over a dominate mode associated with spin fluctuations 
with fl s f = 260 meV (A s /, z = A s /^ = 1.6). Inset: dN/duj 
and c?F(y — An) highlighting the correspondence between 
the low-energy mode and structure in N(lj). (b), (c) A(ui) 
and Z{uj) for this case. The red dashed lines indicate the en- 
ergy of An + fin where fin is the energy of the center of the 
low-energy spectra density (inset, panel (a)). 



where 8<p and SZ are the el-boson contributions to the 
self-energy. The DOS (for u> ^> Ao) can then be written 
as 



N s (lu) 



N 



f 



A 2 
2w 2 



8<j){u) SZ{lu) 
<t>0 Z 



(10) 



If the el-boson contribution to pairing is small, as is the 
case in Fig. 2, 8cf> can be neglected and one sees that the 
fine structure tracks the structure of 5Z(uj) (Fig. 2c). 
In Fig. 2b, this is seen as the suppression of Ai(u) as 
u> — > An + Q and results in a dip structure in N(uj) 
(Fig. 2a) with no pronounced shoulder on the low-energy 
side of the renormalization. We also note that the boson 
energy scale remains as the minima in dN / du (inset of 
Fig. 2a). 

Experimentally, the modulations in the tunneling spec- 
tra appear as a dip-hump structure with no pronounced 
shoulder on the low-energy side of the modulations. ' 22 
We therefore conclude that, within validity of Eliashberg 
theory, the bosonic mode responsible for the LDOS renor- 
malizations cannot be the sole contributor to pairing in 
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FIG. 3: (Color online) The structures in the DOS N(u) and 
gap function A(cj) when the structure of the bare band has 
been retained. (al),(a2) N(tj) and A(ui), respectively, for a 
ci-wave superconductor coupled to a single low-energy mode 
which provides all of the pairing. (bl),(b2) N(uj) and A(u>) 
for the two mode case analogous to the case shown in Figs. 
Icl,c2. The insets of panels (al) and (bl) show dN/dui and 
the low energy component of a 2 F(u + Ao) in order to show 
the correspondence between structure in N(u) and peaks in 
a 2 F(v). (cl), (c2) The real and imaginary parts of Z(u) 
calculated in the infinite- (solid) and finite-band (dashed) 
Eliashberg formalisms. Here, a single iteration of the Eliash- 
berg equations has been assumed with coupling to an Einstein 
mode centered at fi = 52 meV and assuming a ci-wave gap 
with Ao = 35 meV. 



the cuprates otherwise a pronounced shoulder on the 
low-energy side of the renormalization would be present. 
More likely, the low-energy boson renormalizes over the 
dominant source of pairing. This approach has been in- 
voked previously in order to account for the ARPES kinks 
since the value of A z ~ 0.3 — 0.5 needed to reproduce the 
kink is too small to provide enough pairing to account 
for the large gap values. 19,34 

Although the two mode picture can account for the 
lack of shoulder feature in the LDOS in the infinite band 
formalism, the hump structure on the high-energy side 
of the renormalizations is absent. This is due to the 
approximations inherent to the formalism used so far, 
which neglect the role of both bandstructure and the 
narrow bandwidth of the system. In Fig. 3 solutions to 
the Eliashberg equations given by Eq. (4), which retain 
the full k-dependence of the band structure, are shown. 
Again, two cases are considered in analogy to Figs. 1 and 
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2: a single mode with A z = 1.9 = A^, T = 1 meV, and 
il = 52 meV, and a two-mode case with the modes pa- 
rameterized by Slb, s f — 52, 300 meV I\ s y = 1, 30 meV, 
A 2:S / = \<i>, s f = 1.7, A z ,& = 0.52 and A^f, = 0. In order 
to model a realistic bandstructure for the low energy dis- 
persion have assumed a 5-parameter tight-binding model 
obtained from hts to the low-energy dispersion observed 
by ARPES. 44 Fig. 3al, a2, show N{uj) and A(w), respec- 
tively, for the single low-energy mode which is pairing, 
analogous to Fig. 1, while Fig. 3bl,b2 shows the results 
for the two mode calculation, analogous to Fig. 2. 

The first observation is that the shoulder feature re- 
mains when a single mode contributes significantly to 
pairing, while a pronounced dip-hump structure is pro- 
duced when the mode renormalizes over another domi- 
nant mechanism. In the latter case the hump is more 
pronounced due to the additional self-energy contribu- 
tions when scattering to the large DOS near the van 
Hove. In the cuprates the energy of the phonons (and 
spin resonance mode) lie close in energy to both the su- 
perconducting gap Ao and the van Hove energy e(0, Tr/a). 
Because of this near degeneracy of energy scales there is 
an overall enhancement of the self-energy due to the in- 
creased density of states to which the bosons can couple 
anti-nodal region. This degeneracy has also been noted 
in studies on the temperature dependence of the el-ph 
self-energy observed by ARPES. 15,16 To illustrate this 
point, in Fig. 3cl,c2 compares the real and imaginary 
parts of Z(lo) in the two formalisms. (In order to com- 
pare comparable cases Z(uj) as been calculated assuming 
a single iteration of the Eliashberg equations and the re- 
sults have been normalized by the value of A z .) When 
the bandstructure has been retained Z(u>) (as well as 
X and (j)) develops additional structure and Zi(uj) be- 
comes negative at energies larger than E(0, tt)+CIq, where 
E 2 (0, 7r) = e 2 (0, tt) + A 2 (0, tt) is the energy of the quasi- 
particle at the van Hove singularity, resulting in the hump 
structure observed in N(uj). If the bandstructure is ne- 
glected Zi(uj) remains positive and smoothly approaches 
zero as w — > oo. One can therefore conclude that the 
structure in the underlying band can contribute to the 
structure in the self-energy and must be considered in 
realistic treatments of narrow bandwidth systems such 
as the cuprates. 

The second observation made of Fig. 3 is that the 
correspondence between the minima in dN/duj and peaks 
in a 2 F no longer holds for a d-wave superconductor with 
a finite bandwidth (insets of Figs. 3al and bl). For the 
single-mode model the peak in o?F(y) corresponds to 
the minimum in N(lu) or the root in dN/diu while for the 
two-mode model this energy scale does not correspond 
to any feature in dN/duj. For the choice of a 2 F and ek 
used here, the energy scale A + fio is located between 
the root and maximum in dN / du and either energy scale 
gives a reasonable estimate for the mode energy. 

We also note the presence of secondary features in the 
DOS at an energy E(0,ir) + Slo which is absent when 
momentum structure of the band is neglected. The pres- 
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FIG. 4: (Color online) (a) N(uj) and (b) dN/dio for u in the 
neighbourhood of the low-energy boson renormalizations for 
various values of the low-energy mode's spectral width Ft, (in 
meV). 



ence of this feature can complicate the identification of 
the mode energy if the boson spectral density is broad 
enough. In Fig. 4 N(uj) and dN/duj are plotted for 
the two-mode model for various values of the low-energy 
mode's spectral width Tf,. Here, we have kept the details 
of the high-energy mode fixed to the same values used in 
Fig. 3 and adjusted the strength of A& such that the total 
value of A Zj are fixed. As the value of Tf, is increased, the 
width of the structure in N(u) is increased and eventu- 
ally becomes broad enough that the double dip structure 
merges into a single dip-hump feature. As a result, the 
peaks in dN /du merge into a single peak whose maxi- 
mum occurs between Slo + Ao and fio + E(0,ir). The 
maximum in dN / du is very sensitive to the width of the 
boson spectrum as well as the energy of the van Hove and 
is therefore an unreliable indicator of the boson mode 
energy. However, the position of the root in dN/cLu on 
the low energy side of the renormalizations is much less 
sensitive to these factors. Although this feature slightly 
underestimates the mode energy, it does appear to be a 
robust indicator of the mode energy. 

The qualitative change in the phonon fine-structure is 
the central result of this paper. It reconciles the dis- 
crepancy in the ~ 52 meV scale observed in STM 3 and 
the ~ 70 meV "kink" observed by photoemission in Bi- 
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2212. 15 In Ref. 3 the maxima in d 2 I/dV 2 oc dN/dtu was 
taken for the mode estimate, corresponding to the shoul- 
der of the dip-hump structure of the DOS. However, as 
shown above, the energy scale of the mode is mode ac- 
curately given by the minimum (the dip) in N(lu) (or a 
root in d 2 I/dV 2 ), this choice results in an overestimate 
of the mode energy. A closer examination of Fig. 1 of 
Ref. 3 reveals that the minima in the DOS is ~ 15-20 
meV lower in energy. This brings the mode estimate in 
line with the energy of the Bi g modes invoked to explain 
the kink in the nodal region of superconducting Bi-2212. 
Indeed, Ref. 22 has tracked the minimum in the related 
Bi-2223 system, which also has strong coupling to the B\ g 
modes, and obtained a mode energy of 35 meV, consistent 
with our findings. The qualitative difference in structure 
also provides a pathway to experimentally distinguish be- 
tween fine structure due to co-tunneling effects, where a 
maximum in iV(f2o + Ao) oc dl / dV^^rta+ko occurs, and 
intrinsic el-boson coupling where a minimum is expected 
if the mode is not dominating pairing, or a shoulder is 
expected if the mode is contributing heavily to pairing. 
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FIG. 5: (Color online) (a) The density of states calculated for 
coupling to the Bi g phonon branch at T = 0. Each spectra 
is calculated using a 5-parameter tight-binding bandstucture 
and a d-wave gap with Ao = 20 (black), 40 (red) and 60 
(blue) meV. (b)-(d) dN/duj for the indicated gap size. The red 
dashed lines indicate the position of the roots corresponding 
to estimates for Ao and Ao + fio, respectively. 



IV. CONSIDERATIONS FOR BI-2212 

We now turn to a model calculation for Bi-2212. Here 
we consider coupling to the out-of-phase Cu-0 bond 
buckling B\ g phonon branch, previously invoked to ex- 
plain the renormalization in the bandstructure of Bi-2212 
observed by ARPES. 15 Since we are now concerned with 
how the phonons renormalize over a dominant interac- 
tion we consider only a single iteration of the Eliashberg 
equations and treat the phonon mode as a dispersion- 
less Einstein mode with ilo = 36 meV. This treatment 
is identical to that used in previous works examining 
the dispersion kinks observed by ARPES. 15,16 ' 19 For sim- 
plicity, we expand the coupling constant as <?(k, q)| 2 = 
9l +9l Y d{k)Y d (p) and set g z ^ such that A Zj0 = 0.31, 0.1, 
comparable to the values obtained in previous works. 15 
After analytic continuation, the zero-temperature expres- 
sions for the imaginary parts of the self-energies are: 

w z 2 (k, w ) = ^J2\9(K<i)\ 2 s(E P + n -cj) (ii) 
p 

X2(k,w) - -^5>(k !q )| 2 |^(£ p + fio-u,) 
p p 

<fe(k,c) = -^^| 5 (k,q)| 2 ^( J B p + fi - W ) 
p p 

which are evaluated for uj > 0. Here, E(k) = 
•\/e 2 (k) + A 2 (k). The real-part of the self-energies are 
obtained via the Kramers-Kronig relations. In these cal- 
culations the superconducting gap Ao is taken as an in- 
put parameter and we supplement the real part of <f> in or- 
der to maintain the value of the gap at the gap edge. 19,34 
Finally, an intrinsic damping r = 5 meV, independent of 
lo and k, is added to the imaginary part of Z(k.,uj). 



The DOS results are presented in Fig. 5, where we have 
chosen three different values of Ao to mimic the variation 
in the LDOS observed in different regions of an optimally 
doped sample. The calculated DOS for the three gap val- 
ues all show a clear dip in the spectra around Ao + Slo 
as well as the secondary structure at E(0, ir) + Slo due to 
the van Hove singularity. We also note that the features 
associated with the phonon renormalizations are more 
pronounced due to the use of a sharp Einstein mode. Es- 
timates are obtained empirically for the mode energy and 
gap size from the calculated DOS via the gap referencing 
procedure used in Ref. 4. In order to avoid complications 
associated with the van Hove singularity in determining 
the gap magnitude on the occupied side (u> < 0) of the 
spectra we work on the hole side (ui > 0) and identify the 
energy Ao + flo with the root in dN/du. When the gap 
referencing procedure is applied to the data, the result- 
ing gap estimate is equal to the quadrature addition of 
the gap on the Fermi surface Ao and the damping term 
r. The empirically determined mode energy is there- 
fore underestimated since the effective mode position is 
fio + Ao — \/ Aq + r 2 . Since we have used a constant 
value of r for each value of A , the small gap data has a 
larger ratio of r/A , and thus the extracted mode energy 
systematically deviated from fig. 

A constant value of T does not capture the local in- 
homogeneity of the parameters entering into the DOS 
itself. Specifically, we now modify the magnitude of the 
el-ph coupling and inelastic damping T included in both 
the spectral function and evaluation of the self-energy, 
as a function of gap size. We naively associate the larger 
gap with "underdoped" regions which, due to the reduc- 
tion in screening of the el-ph interaction, leads to an in- 
crease in the relative strength of the coupling with gap 
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FIG. 6: (Color online) A waterfall plot of N(tu) calculated 
for doping values spanning the large gap to small gap regions. 
Each DOS was calculated for coupling the same mode used 
in Fig. 5. The gap values indicated in the legend denote 
the input values of Ao. (b),(c),(d) dN/dui for selected DOS 
presented in panel (a) for the hole side (cu > 0) of the spec- 
trum. The red arrows indicate the roots used to estimate the 
energies Ao and Ao + fio • (e) The mode energy estimate ob- 
tained from the position of the local minima relative to the 
coherence peak. The blue data points correspond to data for 
16 O simulations while the red data points correspond to 18 O 
simulations. 



size. At the same time, damping effects are taken to in- 
crease together with the gap size to mimic the crossover 
to smeared gap structures found in the large gap regions. 
This is modeled as an increasing ratio of T/Aq in the 
large gap region. The DOS was then recalculated for 
input gap values ranging from 15 to 50 meV. 

The new DOS spectra obtained are presented in Fig. 
6a. For small gap inputs one can see sharp coherence 
peaks followed by a well defined dip-hump structure as- 
sociated with the el-ph coupling. For larger gaps, the as- 
sociated larger damping smears the coherence peaks and 
they are washed away for Ao = 50 meV. The gap and 
phonon energy scales can now be extracted from dN/duj, 



which are shown for selected spectra in Figs. 6b-d. The 
resulting correlation between the extracted Ao and flo 
are shown Fig. 6e (blue dots). The anti-correlation which 
emerges between the two energies stems from the progres- 
sive underestimation of the mode energy as the gap size 
and r are increased. 

In order to model the isotope effect, these calcula- 
tions were repeated with adjustments appropriate for 
the replacement of 16 with 18 O. Specifically, this in- 
cludes a shift in the phonon frequency by a factor of 
\J Mi 6 /Mi 8 and a decrease in the overall coupling by a 
factor (Mig/Mis) 1 / 4 . (This substitution leaves the val- 
ues of A 2j unchanged.) The red data points of Fig. 6e 
show the correlation between the estimate for f2 and Ao 
obtained for 18 O upon repetition of the previous calcu- 
lations. In both cases, the anti-correlation persists and 
a clear isotope shift can be seen, which is on the or- 
der of that observed and that one would expect based 
on the known shift in the phonon energy. The overall 
agreement with the experimental data is good, and the 
anti-correlation can be accounted for relatively well by 
incorporating damping effects into a simple el-ph picture. 



V. CONCLUSIONS 

We have examined the signatures of el-boson coupling 
in a d-wave superconductor when the boson mode does 
not provide the dominate pairing interaction. The re- 
sulting fine structure for the modulations in the DOS are 
qualitatively different than those found in the conven- 
tional s-wave superconductors. We found that the man- 
ifestation of bosonic energy scales in the LDOS depends 
on an interplay between the strength of the el-boson cou- 
pling in each momentum channel, as well as the details 
of the underlying bandstructure. In order to reproduce 
the spectra observed experimentally by STM we found 
that both the full structure of the band, as well as the 
differing contributions of X z and for the boson mode 
in question have to be accounded for. Taking these con- 
siderations into account, the 52 meV scale reported in 
Ref. 3 is found to be an overestimate of the mode energy 
and instead, the local minima in N(ui) should be con- 
sidered. The mode energy extracted from the minima of 
the spectra reported in Ref. 3 agrees well with the energy 
of the B\g phonon modes as well as a more recent mea- 
surements reported for Bi-22 23. 22 This fact reconciles the 
scales observed by STM with those observed in numer- 
ous ARPES experiments. With the revised estimate for 
the mode energy, we then calculated the LDOS includ- 
ing coupling to the B\ g phonon modes. Using a simple 
consideration for the local intrinsic damping we found 
that the model is able to reproduce both the structure 
of the renormalizations and the observed anti-correlation 
between the extracted values of JIq and Ao . The success 
of the el-ph model in accounting for the data obtained by 
both probes provides further evidence in support of the 
phonon interpretation of the electronic renormalizations 
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observed now in both ARPES and STM experiments. 
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